Supercurrent induced domain wall motion 
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We study the dynamics of a magnetic domain wall, inserted in, or juxtaposed to, a conventional 
superconductor, via the passage of a spin polarized current through a FSF junction. Solving the 
Landau-Lifshitz-Gilbert equation of motion for the magnetic moments we calculate the velocity of 
the domain wall and compare it with the case of a FNF junction. We find that in several regimes 
the domain wall velocity is larger when it is driven by a supercurrent. 



I. INTRODUCTION 

Spintronics has attracted considerable attention both 
from the fundamental poin t of view and from the point 
of view of application^II^. Related to applications in 
magnetic registers, the manipulation of magnetic domain 
walls by spin polarized currents has been a topic of inter- 
est, as opposed to direct manipulation by application of 
local magnetic fields. Domain walls are present in mag- 
netic materials and thin magnetic films. Their function 
is to lower the magnetostatic stray field energy^l The ba- 
sic phenomena of the DW motion occur in a submicron- 
size ferromagnetic strip^. The manipulation of the lo- 
cation of the domain walls results from the action of 
spin torques on the magnetic moments^, due to a polar- 
ized currentPl^. T his ha s been analyzed theoreticall}!^''^ 
and experimentall}ff^^. Typically the magnetic domain 
walls are inserted in magnetic semiconductors but inter- 
est in other types of structures has arisen lately involving 
superconductor^. 

We have been interested in the possible int erplay 
between magnetic moments and superconductors'^^, 
specifically the possibility of ordered magnetic moments 
in the superconductor or in its vicinit}'^^. Magnetism and 
superconductivity typically compete and, therefore, most 
heterostructures considered tend to separate the ferro- 
magnetic and superconducting regions (often separated 
by an insulator to prevent proximity effects). The pres- 
ence of randomly located magnetic impurities in a con- 
ventional superconductor destroys the superconducting 
order for s mall impurity concentrations^^ but we have 
showrPESlthat, if the magnetic moments are correlated, 
superconductivity is much more robust and prevails for 
much higher concentrations. Therefore, if the magnetic 
moments are somewhat diluted we expect that super- 
conducting order should remain. We will consider here 
diluted magnetic moments. 

The effects of the polarization of the electrons due to 
the magnetic moments is well known to be very local. 
Therefore, one may expect that passing a spin polar- 
ized current through the superconductor (as in a normal 
metal) will induce a spin interaction between the mag- 
netic moments and the spins of the conduction electrons. 



In a conventional superconductor, with singlet pairing, 
the spin density dies out inside a clean superconduc- 
tor at a distance of the order of the coherence length. 
Therefore, in the case of an infinite system one does not 
expect a spin torque on the magnetic moments due to 
the incident spin polarized current. However, the size 
of the nanostructure may be comparable to the coher- 
ence length and, therefore, the spin density may not die 
out inside the superconductor. Also, due to the polar- 
ization effect of the local magnetic moments on the spins 
of the conduction electrons, a smaller decay of the spin 
density occurs. We have shown recentl}'^^ that there is 
an appreciable spin torque induced on the magnetic mo- 
ments that is comparable to the one observed in a normal 
metal. Clearly, one expects that considering system sizes 
of the order of the coherence length will enhance the ef- 
fect. Also, spin triplet superconductors will circumvent 
this issue altogether. 

It is therefore interesting to further explore a junction 
of the type ferromagnet-superconductor-ferromagnet 
(FSF) to find regimes where the domain wall motion 
induced by the passage of the supercurrent may be 
more effective in the FSF junction with respect to the 
ferromagnet-normal-ferromagnet (FNF) junction. Here 
we study the motion of the domain wall solving the 
Landau-Lifshitz-Gilbert (LLG) equations for the mag- 
netic moments. We find that often the velocity of motion 
of the domain wall is actually larger for the supercurrent, 
as compared to a normal metal. 

A spin polarized current leads to an accumulation of 
spin density which interacts with the local moments in- 
ducing a torque and consequent rotation. This leads ulti- 
mately to the intended shift of the domain wall position 
by methods involving currents. These are more efficient 
and fast as compared to the application of local magnetic 
fields to flip the magnetic moments. For recent reviews 
see Refs.^^ and references therein. Due to the finite 
size of the systems considered the motion of the domain 
wall is limited by the system size. When the domain wall 
approaches the boundary it gets distorted and stops its 
motion. We are therefore interested here in studying the 
early time regimes during which the domain wall is set 
into motion and has not yet distorted appreciably. 
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II. MODEL 

In this work we consider a junction of either FNF or 
FSF types. For simphcity we consider a one-dimensional 
model system which is a good description of a narrow 
superconducting wire. Inside the superconducting region 
we place local magnetic fields. On the left-hand side the 
ferromagnet exchange field, /i, points in the z-direction 
and on the right hand side the exchange field points in 
the —z direction. We take the domain wall resulting from 
the diluted magnetic moments inserted in the supercon- 
ductor centered at the midpoint of the SC (see Figs. 3,4 
below for t = 0). We use a lattice formulation for a ID 
model system oriented along the x axis, with Hamilto- 
nian: 

H = He + He.s+Hs (1) 

where 

He = - ^ (4aCn+la + H.C.) 
n,(T 

+ E(^"4^ii + H-c-) (2) 

n 

is the electronic part of the Hamiltonian, a =t, | denotes 
the spin projections along the z axis, we set the hopping 
to unity and, thereby set the energy scale, and we choose 
the chemical potential to be zero. 

Hc-s = - E J^„ • s„ (3) 

n,cr,cr' 

is the interaction between the spin density of the con- 
duction electrons and the impurity spins, with Sn = 
^aa'C-l^a^na' whcrc CT = [(7^ ^ CF^ ^ CF^) are the Pauli ma- 
trices, we assume the local magnetic moments Sn to be 
of unit length and to be Zeeman coupled to the electrons 
in the superconductor, behaving as local magnetic fields 
JSn = Jn in the superconductor and 

Hs = -Y.J^-^n-Sn+i + ^Y.(Slf (4) 

n n 

It is convenient to introduce a planar anisotropy ky^ 
which for positive ky favors a state where the spins are 
in the x — z plane (the plane of the initial domain wall) . 
At the superconductor interfaces with the F or N sys- 
tems {Nsl^Nsr) we introduce apotential term, [/, that 
simulates the interface disordei . The local moments 
are distributed evenly inside the superconducting region 
and interact with each other via a nearest-neighbor fer- 
romagnetic interaction J^x- The superconducting region 
includes the magnetic moments that constitute the do- 
main wall, as shown in Fig. 1 of Ref.^^. Often we take 
the number of sites N = 160, and Nsr-Nsl = 100. We 
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FIG. 1: Evolution of < > as a function of J. The various 
plateaus are obtained through quantum phase transitions. 

are considering for simplicity a singlet s-wave supercon- 
ductor. 

The Bogoliubov-de Gennes (BdG) equations^^ for i^c + 
Hc-s determine both the equilibrium electronic states 
and the scattering states when a current is passing 
through. Their solution provide the wave functions in 
the usual way (see for instanc^^. The BdG equations 
are solved self-consistently to find the profile of the gap 
function inside the superconductor, A^. The domain wall 
shape is calculated self-consistently in mean- en- 
suring that the energy is minimized and a vanishingly 
small equilibrium torque. The procedure leads to a do- 
main wall where the magnetization direction interpolates 
between the z-direction on the left hand side and the —z 
direction in the right hand side, according to the orienta- 
tions of the exchange fields in the two ferromagnets. The 
procedure is described in detail in Ref.f^ and in Ref.l^. 

Using the transfer matrix method we determine the 
wave functions at every site in the heterostructure and 
calculate the various physical quantities such as the local 
spin densities, the spin and charge currents and the spin 
torque at each site. The torque at point n is defined by 
the cross product between the magnetic moments and 
the local spin density, s^, 

T(n) = 2JSn X Sn- (5) 

Assuming spin conservation the spin torque may also be 
calculated by the difference in spin currents as 

T^{n)=j^n-l)-j^n), (6) 

where /3 = x^y^z. The expressions for the spin currents 
are given in Ref.^^. 

In experiments, a potential difference, is imposed 
between the two sides of the heterostructure in a stan- 
dard wa}E^H^ leading to a current that moves from left 
to right. Imposing a potential difference at the ends of 
the heterostructure leads to an overall torque on a given 
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FIG. 2: Spatial variation of the gap function for different values of J from J = 0.3, • 
increases by one as the average magnetization jumps from one plateau to the next. 



,2.5. Note that the number of peaks 



magnetic moment, that is the result of an integral over 
the incident energies, up to the applied potentiaP^^. This 
torque leads to the motion of the magnetic moments, 
governed by the Landau-Lifshitz- Gilbert equation^. 



III. EQUILIBRIUM PROPERTIES OF THE 
HETEROSTRUCTURE 



We begin by considering a situation with no current 
flowing through the system. In this case the BdG equa- 
tions yield the equilibrium solution. In Fig. [l]we show 
the average magnetization of the electrons over the left 
half of the superconductor, < >, along the z-direction, 
induced by the coupling to the local spins, as a function of 
J. We consider 15 spins distributed evenly in the super- 
conductor (density 0.15). We see in the superconducting 
case that, as J increases, the magnetization increases, 
but at some points it changes discontinuously between 



various plateaus. These discontinuiti es are due to quan- 
tum phase transitions in the systeni^^^. In the case 
of one impurity and for small J the average magnetiza- 
tion vanishes. Even though the local magnetic moment 
polarizes the spin density of the electrons at the impu- 
rity location, the spin density of the electrons has fluc- 
tuations (like Friedel oscillations) that compensate the 
perturbation. Above a certain critical Jc the system is 
no longer able to shield the perturbation and the over- 
all magnetization jumps discontinuously to 1/2. Various 
other discontinuities occur at this flrst order quantum 
phase transition due to a level crossing, such as the gap 
function changing sign at the impurity location, various 
entanglement^^ measures and the partial state fldelity^. 
Increasing the number of magnetic impurities the discon- 
tinuities are still present, but at very small J the magne- 
tization is not strictly zero^^, as shown in Fig. [l] We also 
show in Fig. [l]the magnetization for the normal case. We 
see that there are no discontinuities and the magnetiza- 
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tion changes smoothly. At smah J the magnetization is 
smaher in the superconducting case. As the first tran- 
sition to a plateau occurs, in the superconducting case, 
the magnetization becomes larger in the superconducting 
case. 

The profile of the gap function is shown in Fig. [2] for 
various spin couplings, J. At small J the gap function is 
approximately constant over the superconducting region. 
It vanishes outside the superconductor and it has very 
small fluctuations due to the finite size of the SC region 
and the finite number of impurity spins. As the coupling 
grows the fluctuations in the gap function increase con- 
siderably and are modulated by oscillations that increase 
in number as J grows. Note that, interestingly, the num- 
ber of peaks in the oscillations of the gap function as a 
function of space along the chain is related to the various 
plateaus of the total magnetization. As one goes from 
one plateau to the next the number of peaks increases 
by one. This is reminiscent of a LOFF^^ state where the 
gap function has real space oscillations due to the finite 
momentum of the Cooper pairs resulting from the split- 
ting of the up and down spin bands due to a magnetic 
field. We recall that in a superconductor with magnetic 
impurities, if J is large enough, the gap function changes 
sign at the impurity locations. Here, however, we have an 
heterostructure. Interestingly, we find that even though 
the gap function decreases significantly at the impurity 
locations as J grows, it only changes sign at very large 
J ~ 2.5. This is in contrast with an infinite superconduc- 
tor where the gap function changes sign at the location 
of the first quantum phase transition. This is shown in 

Fig.m 



IV. LANDAU-LIFSHITZ-GILBERT EQUATION 
AND DOMAIN WALL MOTION 

Consider now that the spin current is turned on. The 
LLC equation for the motion of a magnetic moment S in 
the presence of an effective magnetic field H and subject 
to a spin current (Slonczewski term^) giving origin to 
a spin torque r is written as 



~dt 



= --fS X H 



a dS 



(7) 



Here S is the length of the magnetic moment, 7 is the 
gyromagnetic ratio and a is the damping constant. The 
magnetic field ii" is a generic notation for an effective 
field defined as a derivative of the energy with respect 
to the magnetization (spin). Due to the spin torque the 
magnetic moments in general evolve out of the plane and 
therefore we assumed the system has an anisotropy term 
that favors the plane. Also, the interaction between the 
local magnetic moments leads to a term that contributes 
to the effective magnetic field. The effective magnetic 
field can be written as 



The contributions of these two terms to the effective mag- 
netic field are 



and 



— kySyy 



Hf{n) = Je^Sp{n + l) 



(9) 



(10) 



where f3 = x^y. 

Taking the cross product of the LLC equation with 
the magnetic moment and using that its length is fixed 
{S • dS/dt = 0), the LLC can be reformulated as 



dt 



1 



S 

a 1 



1 + a2 1 + ^2 5 



S XT 



where 



7L 
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(11) 



(12) 



In terms of the spherical angles 0^ cp it is written as 
dO ( ~ ~ \ 1 1 la 



1 a 



T9 



(13) 



where Hq^ tq and H^^ are the spherical components 
of the vectors H and r, respectively. The anisotropy 
and the exchange term imply extra terms in the LLC 
equations of the form 



= —kyS sin cos sin (p 
= — kyS sin^ simp cos (p 



(14) 



and 

rrex 

- 

rrex 



Jex [cos9„ sin6l„+i cos {ip„ - (pn+i) - sin^„cos6l„+i] 

-Jex[smen+lSm{ipn - fn+l)] (15) 



H = Hy + w 



(8) 



These equations are solved iteratively: for a given do- 
main wall configuration, at a certain time, we calculate 
the electronic properties, such as the spin torque from 
Eqs. |5|6[ This torque is then used to evolve the magnetic 
moments to the following time using the LLC equation. 
The LLC equations change the orientation of the spins. 
This new configuration is then used to calculate the new 
spin torque and so on. To solve the LLC equations we 
used a standard second order Runge-Kutta method or 
the Heun method, yielding similar results. 



V. EVOLUTION OF THE DOMAIN WALL 

We start with a domain wall that interpolates between 
opposite exchange fields in the ferromagnetic regions. 
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t=0. t=0.5 




FIG. 3: Time evolution of a domain wall of 15 spins with parameters J = 1, Jex — 0.5, V — 1. At the initial time (t = 0) the 
domain wall magnetic moments are contained in the x — z plane. As time evolves the location of the center of the DW shifts 
to the right. 



t=0. t=0.4 




FIG. 4: Time evolution of a domain wall of 41 spins with parameters J = 1, Jex — 0.5, V = 1. 



6 




FIG. 5: Time dependence of the (a) distance between the center of the DW, xdw, and the midpoint of the SC and (b) velocity 
of the DW for different values of J from J = 0.1,---,2.5. In black (full lines) are the results for the superconductor and in red 
(dashed lines) the results for the normal system. 
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FIG. 6: Dependence of the domain wall velocity with the various parameters. We consider the dependence with (a) the 
interface disorder, [/, (b) the exchange between the spins, Jex and (c) the applied potential, V . In panel (d) and we compare 
the velocity for the superconductor and the normal metal for the parameters J — 2, Jex = 1.5, F = 0.5, a = 0.5. 



Initially the magnetic moments form a Neel-type domain 
wall in the x — z plane. As time evolves the domain 
wall moves to the right due to the spin torque exerted by 
the spin polarized current. The spin torque has in-plane 
and out of plane components that rotate the magnetic 
moments. The combined effect of the spin torque, the 
exchange coupling between the magnetic moments, the 
magnetic anisotropy and the damping term determines 
the time evolution of the spins. The magnetic moments 
tend to rotate to align with the spins on the left-hand 
side (leading to the intended displacement of the center 
of the DW) but also tend to move out of the x — z plane 
(this effect is counteracted by the magnetic anisotropy 
and by the damping term). 

In Figs. [3] and |4] we show the time evolution of the 
domain wall for different parameters. There are several 
parameters that need to be specified which implies a very 
large parameter space. As a standard case shown in Fig. 
[3] we take J I, Jex = 0.5, V = l,ky = 4,a = 0.02, /7 = 
0,7 = 2.2, = 0.2. The parameter h is the exchange 
magnetic field in the ferromagnetic regions. We take h — 
0.2 in the left F and = -0.2 in the right F. 

Here we are interested in studying the early time evo- 
lution, up to the moment the domain wall has come close 
to the edge of the spin configuration. In particular, we 



are interested in comparing the time evolution using a 
supercurrent with the one obtained from a standard nor- 
mal metal. In Fig. |5^ we show the displacement of the 
center of the DW as a function of time for various pa- 
rameters. We define xbw as the distance between the 
center of the DW and the midpoint of the SC. 

To better quantify the motion of the domain wall we 
calculate the velocity shown in Fig. [5]3. In the initial 
configuration the spins on the left hand side have an ori- 
entation such that 6> < 7r/2 while on the right hand side 
^ > 7r/2. As the domain wall moves to the right, the ori- 
entation of spins on the right hand side will successively 
cross 7r/2. We define the center of the DW to be the spin 
with Q = 7r/2. The velocity of the DW is {dxDw/dt).. 
This is shown in Fig. ^jp for different values of J, where 
we compare the velocities as a function of time for the 
superconductor and a normal system. This last case is 
simply obtained turning off the superconducting pairing. 

Typically, at early times, the velocity is small. It in- 
creases over time and then starts to saturate. This is 
related to either basically all spins have already flipped 
or the DW starts to get distorted and the last spins do 
not flip at all. As discussed in Ref.Hlthe spin torque has 
in general a component that leads to the motion of the 
magnetic moments out of the original DW plane. This 
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t=1.25 t=1.5 




FIG. 7: Time evolution of domain wall of 15 spins and J = 1, J ex — 0.5, = 1 for later times. For the parameters considered 
here the DW center point reaches the right hand side border and the remaining non-aligned magnetic moments move out of 
the X — z plane. The boundary condition imposed by the right hand side ferromagnet fixes the orientation of the last magnetic 
moment (t = 1.25). At later times t = 1.5, 1.75 the DW starts to distort and several spins rotate out of the original plane. 



t=0. t=0. 75 




FIG. 8: Time evolution of domain wall of 15 spins and J = 2, J ex — 1.5, y = 0.5, a = 0.5. Increasing the couplings between 
the magnetic moments and their coupling to the electron spin density, and the Gilbert damping, the domain wall distortion is 
smaller, even though some spins start to move out of the initial plane when the center of the DW approaches the right hand 
boundary (t = 1.65). 
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term implies an acceleration of the motion of the domain 
wall that leads to an increased velocity. The velocity then 
saturates both due to the presence of friction originated 
by the damping term, and by shear finite size effects due 
to the approach to the boundary of the system. As a 
function of the coupling between the magnetic moments 
and the electron spin density, J, we see various regimes. 
At small J the DW moves faster for the superconductor, 
between J = 0.7 — 1.2 the normal metal is more efficient, 
and then it changes from one to the other. It may seem 
therefore that the velocity may be larger if the magneti- 
zation is smaller. However, there are regimes for larger 
J for which the velocity is larger in the superconducting 
phase, such as J = 1.5,2. Also, one would expect that 
the velocity should be larger if the spin torque is larger. 
The values of the velocities are similar for all cases but 
for small J and for the case of J = 1.5 the velocity in the 
superconductor is particularly larger with respect to the 
normal metal. For this particular value of the coupling 
and for this system size the velocity does not seem to 
saturate before the domain wall gets distorted. 

To establish in greater detail the influence of the vari- 
ous parameters we consider a few cases. Due to the large 
parameter space we fix the various parameters at the 
standard case and change one of the parameters to see 
its influence. In Fig. [6] we change the interface disorder 
f/, the coupling between the magnetic moments, Jex^ and 
the external potential, V. As one expects, the velocity is 
typically larger for a clean interface, or at least for small 
interface disorder. The velocity increases as the exchange 
Jex increases and as the potential, V, increases. Small 
interface disorder is important to have larger currents, 
and larger potential has the same effect. Larger currents 
lead to larger spin torques. The effect of the exchange 
interaction is also clear because it tends to collectively 
move the spins. 

From experimental results one knows that the do- 
main walls either stop, distort or change into vortex 
configuration^^^^, and the process has to be restarted. 
This is not related however to the finite size effects of the 
wire. In any case we were mostly interested in describing 
the early times displacement of the DW. 

Since the system sizes considered are finite, the motion 
of the domain wall is conditioned by the end of the sys- 
tem. Typically the domain wall moves from left to right 
and, when the spins close to the edge move appreciably, 
there is an important distortion of the domain wall since 
the magnetic moments rotate out of the original plane. 
This occurs near the right edge of the system but of- 
ten it propagates to the left. In some cases it destroys 
the domain wall on the right hand side and the spins on 
the left hand side may start to precess and move out of 
the plane. Clearly, these situations are not what is in- 
tended since we want to move the center of the domain 
wall but maintain the spins aligned in the original plane 
(defined as x — z). This is illustrated in Fig. [7| Beyond 
the point for which the center of the DW approaches the 



right boundary the spins on the right hand side start to 
disorder (in the sense they loose the alignment with the 
neighboring spins) and the spins on the left hand side 
start to move out of the plane. 

This problem may be avoided (for not very long times) 
by changing the parameters. We present in Fig. [8] the 
time evolution of the same domain wall with 15 spins (re- 
call that the system size is 100 sites and therefore 15 spins 
is rather dilute) where we have changed some of the pa- 
rameters to a set with J = 2, J^^ = 1.5, V = 0.5, a = 0.5. 
The motion of the domain wall is now much better con- 
trolled, at least up to a time of the order oft = 2. Clearly, 
when all the spins have been flipped one should turn the 
current off. We have found that for these parameters it 
is important to increase the coupling between the mag- 
netic moments and the spins of the conduction electrons, 
to increase the coupling between the magnetic moments 
and to increase the damping term, so that the motion of 
the magnetic moments does not lead to a large distortion 
of the DW. As shown in Fig. [6]i in this case the SC veloc- 
ity is considerably larger than the velocity in the normal 
phase. 

Another way to make the evolution of the domain wall 
well controlled is to increase the density of the magnetic 
moments as shown in Fig. [4j The evolution is slower but 
well controlled up to t = 2, for which the domain wall 
has not yet reached the border. 



VI. CONCLUSIONS 

Following previous works where we considered the pos- 
sibility of the effects of magnetic impurities, immersed or 
in the vicinity of a conventional superconductor, paying 
attention to both the effects of the impurities on the su- 
perconductor and of the superconductor on the magnetic 
impurities, we have solved the LLC equations for the 
motion of a domain wall due to the passage of a spin 
polarized current through a heterostructure of the type 
FSF. The new aspect is that the magnetic moments are 
distributed in the S region. Consistently with previous 
results, where the spin torque exerted on the magnetic 
moments may be larger in the superconductor as com- 
pared to a normal metal, we have found that the velocity 
of the motion of the domain wall in the superconduct- 
ing phase may be enhanced with respect to the normal 
phase. 
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